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Abstract 

The main goal of this paper is to provide a fractional stochastic differential équa- 
tion modelling the physical phenomena governed by the Langevin équation in 1- 
dimension. A generalized équation leaning on the fractional Brownian motion (fBm) 
will be proposed, the later will allow a description of the complexity of the physi- 
cal Systems which escape any prédiction of the of the standard Langevin équation. 
We shall begin at first to remind the basic notions of the standard Brownian mo- 
tion (Bm) and the fractional Brownian motion (fBm), then, we shall establish a 
generalization to long memory of the Langevin équation. 

Résumé. 

L'objectif de ce travail est de donner une équation différentielle stochastique frac- 
tionnaire modélisant les phénomènes physiques gouvernés par l'équation de Langevin 
en dimension 1. Une équation généralisée s'appuyant sur le mouvement brownien 
fractionnaire mBf sera proposée, elle permettra une description de la complexité 
des systèmes physiques qui échappent à toute prédiction de l'équation de Langevin 
standard. On commencera dans un premier temps par rappeler les notions de base 
du mouvement brownien standard mB et le mouvement brownien fractionnaire mBf, 
puis, nous établirons une généralisation à mémoire longue de l'équation de Langevin. 
Ce travail s'inscrit dans le cadre de la continuité de l'application du travail de Lounis 
(2007) inspiré des travaux de Nualart (2003). 



1 La mémoire courte et longue 



Nous connaissons de nombreux modèles stochastiques à mémoire courte: on peut 
citer les variables indépendantes, les variables m-indépendantes, certain processus 
de Markov, les moyennes mobiles, la plupart des processus autorégressifs à moyenne 
mobile (ARMA) et de nombreux processus linéaires. L'avantage de la courte mémoire 
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est qu'elle induit, le plus souvent, de nombreux théorèmes limites comme les lois des 
grands nombres, des théorèmes centraux fonctionnels et des théorèmes de grandes 
déviations. Lorsque la mémoire est longue, la situation est souvent différente. Il y a 
au moins un siècle que les astronomes ont noté l'existence de séries d'observations à 
mémoire persistante. Depuis lors, des phénomènes de même type ont été observés, 
notamment, en chimie, en hydrologie, en climatologie, et en économie. De telles 
séries posent évidemment d'intéressants problèmes statistiques et leur modélisation 
ainsi que leur tr aitement stat istique ont déjà fait l'objet de plusieurs travaux de 
recherches. Dans iBeranI (jl994l ). on trouve quelques exemples historiques traités par 
modélisation brownienne fractionnaire. La modélisation de certains phénomènes 
physiques, dits à mémoire longue, peut s'effectuer par l'introduction de termes 
intégro-différentiels à noyaux faiblement singuliers (càd localement intégrables, mais 
pas nécessairement continus) dans les équations de Hamilon - Jacabi ou Euler - La- 
grange qui régissent la dynamique physique. Il existe également des applications à 
la modélisation en chimie des polymè res ou à l a mod élisation de la dynamique à 
l'interface de structures fractales: voir iMehautd (120001 ) pour l'aspect physique ap- 
pliquée et, par exemple, iGiona and Romani (Il992l ) pour l'aspect physique théorique. 
De plus, la dérivation fractionnaire peut paraître naturellement lorsqu'un phénomène 
dynamique est fortement conditionné par la géom étrie du problè me: un exemple sim- 
ple, très instructif, est présenté dans la référence iTorvikl (11984 ). 



2 Préliminaires 



Dans toute cette partie, {Q, JF, P) est un espace de probabilité, {J-'t)t^o une filtration 
sur cette espace, {Xt,t ^ 0} un processus adapté à valeurs dans R^, {d ^ 1), 
{Bt,t ^ 0} est J^t mouvement brownien. [0,T] est un intervalle borné dans M sur 
lequel on définit une subdivision (tj) où. i & {1, . . . , n} et to = < ti < ■ ■ ■ < tn = T , 
et on pose A = max(tj — tj_i). 



2.1 Mouvement brownien ordinaire 



Un mouvement brownien ordinaire B = {Bf, t ^ 0} (ou processus de Wiener ) est un 
processus sur {fl, T , Tt, P) adapté, continu et admettant les propriétés suivantes : 

- Bs) est une variable réelle, de loi gaussienne, 
Bs) est indépendant de JF„ oii m ^ s. 



• p(5o = 0) = 1 P-p.s sur Q. 

• Vs ^ et Vt ^ avec t ^ s; {Bt - 
centrée et de variance t — s. 

• Vs ^ et Vt ^ avec t ^ s; (Bt - 
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2.2 Quelques propriétés du mouvement brownien standard 



2.2.1 Variation quadratique du mouvement brownien standard 

Définition 2.1 La variation quadratique du mouvement brownien 

B = {Bt,t ^ 0} est la limite dans de J2i'=i \Bt^ — i^t. jp quand A tend vers 

0. 

Théorème 2.1 La limite ci- dessus existe dans et vaut T. 



2.2.2 Propriété de Markov 

Définition 2.2 La fonction aléatoire X = {Xt,t ^ 0} est un processus de Markov 
par rapport à T , si pour s ^ 0, t ^ et tout borélien A de on a: 

P {Xt+s G A\rt) = P {Xt+s G A\Xt) . (2.1) 



Cette dernière égalité s'appelle propriété de Markov, cela signifie que la loi de Xt+s 
sachant Tt ne dépend que de Xt. "La loi du futur ne dépend du passé que par le 
présent" . 

Remar k 2.1 Pour plus de déta i ls sur les démonstrations, on peut se rapporter au 
livre de ÏT. Meure et F. Comets I ImM). 



3 Mouvement brownien fractionnaire 



Un mouvement brownien fractionnaire = {Bj^ ,t ^ 0}, de paramètre H dans 
(0, 1), est un processus gaussien centré vérifiant les conditions suivantes : 



(1) B^ est un processus stationnaire; 

(2) E(5f ) = t^^; 

(3) = 0. 

Le paramètre H est l'indice de Hurst. 



Remark 3.1 (1) Si H = ^, B^ est le mouvement brownien standard. 
(2) Si H ^ \, les accroissements du mBf ne sont pas indépendants. 
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3.1 Quelques propriétés du mBf 



Propriété 3.1 B admet comme fonction de covariance, la fonction Rh définie 
pour tout (s, t) G x ]R+ par: 

RH{t,s) = \{e" + s'" -\t-sr). (3.1) 

Remark 3.2 Si H = ^, alors: 



EiBfB^) = Ut + s-\t-s\) 



2 

s 

sAt 

CoviBt, Bs 



On retrouve ainsi la covariance du movement brownien standard. 

Propriété 3.2 La fonction de covariance Rh définie de x M+ dans M est con- 
tinue, semi-définie positive et symétrique. 



3.2 Corrélations à long terme 



Le mouvement brownien fractionnaire est caractérisé par la présence de corrélations 
à long terme lorsque H est différent de |. Un exposant supérieur à | révèle des 
phénomènes de persistance, c'est-à-dire que l'évolution du processus tend à suivre 
ses tendances. Si le processus a augmenté précédemment, la probabilité est forte 
qu'il continue à le faire. Les processus persistants ont une mémoire à long terme, 
c'est-à-dire qu'il existe une corrélation à long terme entre les événements actuels 
et les événements futurs. Chaque observation porte la mémoire des événements qui 
l'ont précédé. A l'inverse, un exposant inférieur à ^ révèle un phénomène d'anti- 
persistance, dans ce cas les accroissements successifs tendent à être négativement 
corrélés. Une augmentation de la variable tend à être suivi d'une diminution, et 
inversement. 



4 Représentation d'un mBf sur un intervalle 



Dans cette partie, la représentation du mBf sera faite sur un intervalle [0,T], une 
technique qui s'appuie sur les noyaux de carré integrable. Fixons l' intervall e [0, T ] 
et considérons les deux cas: H > ^ et H < ^. Premier cas: H > ^ iNualartl (120031 ) 
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avait prouvé que: 

Bf= f KH{t,s)dB, 
Jo 

Où la fonction de carré intégrable Kh définit par: 

Knit, s) = chs^-^ f\u - s)^-iu"-^du, 

J s 



avec t > s, et en clioisissant ch 



H{2H-1) 



, on a: 



rtAs 

/ KH{t,u)KH(s,u)du = RhU, s). 
Jo 



Notons que: 



{t-s 



Deuxième cas: H <^ iNualartI fl2003h avait prouve que: 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



Jo 



Oii la fonction de carré intégrable Kh définit par: 



KH{t,s) = Ch 
avec t > s, et en choisissant ch 



2H 



{l-2H)l3{l-2H,H+^) 



on a: 



rtAs 

/ KH{t,u)KH{s,u)du = Rnit, s) 
Jo 



Notons que: 



^(M) = cMi7-2^^. 



{t - sr 



(4.5) 



(4.6) 



(4.7) 



(4. 



5 Principe d'invariance 



Soit ^ 1} une suite de variables aléatoires réelles indépendantes et de même 
loi, avec = et E^^ = cr^. On considère la marche aléatoire: 



5'n = ^1 H h ^„ , n ^ 1. 



(5.1) 



Le théorème de la limite centrale implique que converge en loi vers J\f{0, 1) quand 
n oo. Derrière cette convergence de variables aléatoires se cache une convergence 
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de fonctions aléatoires. Définissons la ligne polygonale X" extrapolant la marche 
aléatoire S: 

+ , teR^ (5.2) 

X" est une fonction aléatoire réelle continue. On calcule facilement ses accroisse- 
ments: 



loi7V(0,ti) (^7V(0,t2 -ti) 

si n — > oo, par la théorème de la limite centrale. Plus généralement, pour une 
subdivision arbitraire = < ^2 < • • • < on a: 

jim [X'^(t,))r=i - = [B{U) - B{U.,)]T=i (5-3) 

en loi, et donc la convergence en loi des martingales fini-dimensionelles: 

\im[X-{U)]Z^ = [B{U)]Z,. (5.4) 

On admet le théorème suivant: 

Théorème 5.1 (Donsker (1951): Principe d'invariance) La suite de proces- 
sus aléatoires {X"',n ^ 1} converge en loi vers le mouvement brownien. En d'autres 
termes, soit Pn la loi de X'^ sur C(M+,M), et Pb la loi de B. Alors, 

Pn Pb étroitement 



6 Équation de Langevin 

On considère le modèle d'une particule de masse m qui se déplace sur une droite, 
soumise à une force de friction et une force dynamique extérieure Ft qui modélise les 

chocs avec les autres particules. On suppose que la particule n'est soumise à aucune 
autre force extérieure. D'après la seconde loi de la dynamique non relativiste de 
Newton dans un référentiel Galillien, la position xt est alors solution de l'équation 
différentielle: 

mx'l^-hx[ + Ft, (6.1) 
qu'on peut écrire aussi sous la forme: 

mv[ = -bvt + Ff (6.2) 
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avec 6 > 0, le coefficient de friction et Xt la position instantanée de la particule. Si la 
force extérieure est dû à des chocs, nombreux et petits, le principe d'invariance de 
la section précédente nous incite à considérer que /o FJs = aB t, pour un coefficient 
0" > 0, paramétrisant l'intensité de chocs et leurs amplitudes. La force F est alors un 
bruit blanc i.e. la jj dérivée n d'un mouvement brownien, et l'équation ci-dessus est 
improprement écrite, elle doit porter sur les différentielles stochastiques. La vitesse 
Vf = x[ est solution de l'équation différentielle stochastique: 

mdVt = -bVtdt + adBt (6.3) 

au sens où, presque sûrement: 

Vt-Vo+ r —V{s)ds = —Et , Vt ^ 0. (6.4) 
Jo m m 

Cette équation différentielle stochastique fut proposé par Paul Langevin en 1908. 
6.1 Équation de Langevin en dimension 1 

Soit B un {J-'t)t^o mouvement brownien unidimensionnel. On admet la proposition 
suivante: 

Proposition 6.1 La solution de l'équation de Langevin 3\) . partant de Vq G 
L^(jFo), est donnée par le processus d'Ornstein-Uhlenbeck, défini par: 

Vt = + fe-'^'-'^-dB,. (6.5) 

JO m 

On peut montrer par un simple calcul que: 

EVt = e-™*E1/o- (6.6) 



6.2 mBf et l'équation de Langevin 



Si on examine de prêt l'équation de Langevin, on est parti d'une force aléatoire 
Ft qui est approximé par un bruit blanc. Vu que le mouvement brownien standard 
est un processus de markov (par rapport à sa filtration), et pour mieux mettre en 
évidence l'influence du passé, il sera judicieux d'élaborer un modèle qui tient compte 
de l'effet du passé tout en gardant les caractères gaussiens et la stationnarité. A cet 
effet, on intro duit le mouvem ent brownien fractionnaire (processus non markovien 
pour H IStroockl (jl982l ) avait montré qu'on peut approcher un bruit blanc par 

un processus {9e)e>Q- Quand e — cette convergence est en loi dans Co([0,T]). 
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Si > |, alors Jq K^it, s)9^{s)ds converge en loi quand e > vers . 

Si H < i, on a le même résultat de convergence en loi mais avec des conditions 

supplémentaires qu'on notera dans toute la suite (CS). Pour plus de détail s ur les 



conditions de type (CS), on se réfère aux travaux de iBardina et al\ ( l2003l ). Par 



exemple, quand on approche le bruit blanc par le processus: 



s 



ee{s) = Y.^J[k-l,k){^), (6.7) 



k=i 



oii {Ç,k,k ^ 1} est une suite de variables aléatoires réelles indépendantes et de 
même loi, avec E^i = 0, Var^^ = 1, on suppose que les variables aléatoires ont 
des moments d'ordre m et m > j^. 

Dans le but d'écrire l'équation de Langevin généralisée, on traitera séparément les 
cas oii H est supérieur à | et if est inférieur à |. 
Premier cas: H > ^ Considérons l'équation de Langevin (16.30 : 



-bvt + Ff 



(6.8) 



Le processus 9f(t), où OAt) est u n bruit blanc quand e — > 0, (l'existence de ce 
processus est dû à (IStroockl (119821 )). On suppose que pour tout e > 0, on a: 



(6.9) 



oii est une fonction déterministe telle que g,, ^ quand e — 0. Multiplions les 
deux membres de l'équation (16. 3p par la fonction de carré integrable Kh définie par: 



(6.10) 



avec t > s et ch 
pour tout e > 



H(2H-1) 



l3{2-2H,H-\) 



. En intégrant sur l'intervalle [0,t], on obtient 



m / Knit, s)dvs = —h / Knit, s)Vsds + / Knit, s)9^{s)ds — / Knit, s)ge{s)ds 
Jo Jo Jo Jo 

(6.11) 

Au passa ge à la limite e ~^ ^, fn (t, s)9 As)ds converge en loi vers B^^ dans 
Co([0,T]) iBardina et al\ (120031 ): iTudoii (120031 ). Dans ce cas, presque sûrement on a 
l'équation suivante 



m KH{t,s)dVs = -h Knit, s)Vsds + aBf 
Jo Jo 



(6.12) 



L'équation (I6.12p sera utile pour définir l'équation de Langevin généralisée. 
Deuxième cas: H < \ D'une manière analogue en utilisant le noyau 
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Knit, s) = ch 



Z Js 



(6.13) 



avec t > s et en choisissant ch 



2H 



{l-2H)l3{l-2H,H+\) 



type (CS), on obtient une équation analogue à (16.121) . 



, et sous les conditions de 



7 Définition de la vitesse fractionnaire 



Premier cas: H > ^ Après approximation du terme /q Kj^^t, s)Vsds par une somme 



(au sens de Riemann) et vu que KHit,s)ds admet la même unité en temps que 
t^~2^ il est naturel de définir une vitesse fractionnaire en se basant sur l'intégrale 
Jq Knit, s)Vsds normalisée par une fonction de type ^nit) = AHt2~^, où Ah est 
une constante qui dépend de H, dont le rôle sera explicité plus loin. On définit ainsi 
la vitesse fractionnaire V/^ par: 



V,^ = Vo + ^H{t) f Knit, s)Vsds 
Jo 



(7.1) 



Posons fuit) = Knit, s)Vsds et calculons sa dérivée par rapport à t: 



dt 



lim 



\( f KH{t,s)V.4s- f \Hit-h,s)VJs 

Il \J0 Jo j 



h 

dt 



■Vgds + lim 



1 f* 



t~h 



Knit, s)Vsds 



D'autre part 



lim (l r KH{t,s)Vsds] =0 



(7.2) 



On en déduit 



dfH{t) 
dt 



dKnit, s] 
dt 



(7.3) 



Par suite 



dt 



dt 



(7.4) 
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Soit donc 



'dKH{t,s) 1 djVf - Vq) <l>'^{t) 



/o dt dt 



dt 



Multiplions les deux termes l'équation ( 16.12^ par on obtient ainsi 



m$H(t) f KH{t,s)dVs = -b<t>H{t) f KH{t,s)Vsds + a^H{t)Bf 
Jo Jo 



soit 



Jo as 



On introduit la notation suivante: 
Notation: 

Pour toute fonction u{x) = /q ^^-Qi^g{y)dy, on note: 



P^(u(x)) = £f{x,y)^dy 



(Le rôle de consiste à permuter les dérivées partielles). 
Puisque 



dt Jo dt 



VJs 



il est claire que: 
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Finalement, l'équation de Langevin, quand H > ^, s'écrit sous la forme: 
Deuxième cas: H < ^ 

Par un raisonnement analogue, en considérant le noyaux de carré integrable Kh 
correspondant et sous les conditions de type (CS), on obtient l'équation (17.101) . En 

résume, l'équation généralisée de Langevin sur l'intervalle [0,T] s'écrit: Pour H ^ ^ 

Pour H = ^, l'équation de Langevin est donnée par: 

mdVt = -bVtdt + adBt (7.12) 



Remark 7.1 En s 'appuyant sur l'équation ^6.6\) . et en choisissant les noyaux cor- 
respondants avec les conditions correspondantes, on a: 

EV/" = E{Vo){l + <!>H{t) f KH{t,s)e-^'ds) (7.13) 



8 Indice de Hurst et l'analyse R/S 



L'indi ce de Hurst peut être calculer en utilisant la méthode de l'analyse R/S [Hurst 



( 1l95ll ). Pour une série temporelle, X = Xi,X2, . . ■ Xn, la méthode R/S est comme 



suit: 

(1) On calcule la valeur moyenne m: 

1 " 

m = -J2X, (8.1) 

(2) On calcule la série Y de l'ajustement moyen: 

Yt = Xt-m, t = l,2,...,n (8.2) 

(3) On calcule la série Z de la déviation cumulative: 

t 

Zt = Y.Y„ t = l,2,...,n} (8.3) 
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(4) On calcule la série des étendues R: 



Rt = max(Zi, Z2,..., Zt) - min(Zi, Z2, . . . , Zt), t = 1,2, . . . ,n 
(5) On calcule la série R de la déviation standard: 



S, 



\ i=i 



1 * 

Isa-.. 

1=1 



1,2,. ..,n 



(6) On calcule la série R/S des étendues normalisées: 

Rt 



{R/S)t = ^, t = l,2,. 



n 



Hurst iHurstI (Il95ll ) avait montré que lorsque t est assez grand, on a: 



iR/S)t = Xt 



H 



•4) 



15) 



(8.6) 



(8.7) 
(8.8) 



Soit donc: 

ln{{R/S)t) = Hln{t) + ln(A)oùA > 

En faisant une régression simple 011 la variable à expliquer est ln{{R/ S)t) et la 
variable explicative est ln(t ), la pe nte de la droite de régression fournit une estimation 
de l'exposant H de HurstI ( 1951 ). Bien évidemment, Dans la littérature, d' autres 



forme s de cette méthodes existent et pour cela on peut se référer au tr avail d e iHurst 



( )l95ll ). Une estimation de la moyenne de l'exposant H a été donnée par iFellerl (jl95ll ). 



9 Estimation de Ah 



Reprenons l'expression de la vitesse fractionnaire V/^ quand if 7^ |. 

V,^ = Vo + ^H{t) f KH{t,s)Vsds, (9.1) 



ou 



On a donc: 



$^(t) = AHé-"". (9.2) 



J^KH{t,s)VJs ^'-'^ 

Après estimation de H, on utilise le noyau correspondant Kh, puis on estime Ah, 
de la manière suivante: 

On définit sur l'intervalle [0,t], une subdivision < ti < t2 < ■ ■ ■ < = ^- On 
choisira to = et ti très proche de 0. Dans ce cas on estime Ah par: 
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A 



H — — 

n 



1 



n 

E 

1=1 



tr-^(\/,f-K.) 

S',^KH{U,s)Vsds 



Où V£: est la vitesse fractionnaire estimée à l'instant tj. 

Le rôle de Au consiste à mieux approcher les valeurs expérimentales des valeurs 
théoriques. 



10 Utilisation du modèle 

On considère la subdivision précédente de l'intervalle [0,t], et on mesure expérimenta- 
lement les vitesses fractionnaires V^f , à tout instant (tj)^^]^. En utilisant la méthode 
R/S, on estime le paramètre de Hurst H. Avec la valeur de H estimée, et en con- 
sidérant le noyau correspondant, on estime Ah- En revenons à l'équation (17. Il) , et 
en remplaçant H et Ah par leurs valeurs estimées, on obtient une expression de la 
vitesse du phénomène physique considéré. Bien évidemment, la valeur de if, car- 
actérise la série chaotique correspondant au phénomène physique étudié. 



11 Simulation de la vitesse fractionnaire 

Une simulation a été faite pour certaines valeurs de H quand H > \ et H < ^ pour 
une valeur fixée de Ah voir figures . 



12 Conclusion 

Les mouvements browniens fractionnaires (mBf ) modélisent de très nombreux phéno- 
mènes auto-similaires, ou du moins localement auto-similaires: fluctuations boursières 
en mathématiques financières, trafic sur Internet, textures en analyse d'images, ap- 
proches de la turbulence en physique des fluides et la hste est longue. Un domaine en 
pleine expansion est celui des équations différentielles stochastiques régies ou guidées 
par un mouvement brownien fractionnaire. Les expérimentateurs cherchent souvent 
à disposer de modèles leur fournissant des ensembles de fonctions génériques qu'ils 
pourront mieux comparer à leurs observations qu'un exemple particulier de fonction. 
Toutes ces raisons expliquent sans doute la popularité des modèles oii l'irrégularité se 
conjugue avec l'aléatoire. On fait alors appel au probabiliste pour construire des pro- 
cessus (fonctions dépendant du hasard) tels que l'on puisse spécifier leur irrégularité 
sans pour autant pouvoir prédire la forme particulière de la fonction générée. Pour 



13 



illustrer cette manière de penser, et connaissant les propriétés du mBf, il apparaît 
comme un bon processus pour réaliser une telle modélisation. 



H=0.55, A^=1/30 
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Fig. 1. Simulation des vitesses standard et fractionnaires pour H > ^. 
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Fig. 2. Simulation des vitesses standard et fractionnaires pour H < ^. 
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